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ABSTRACT 

"Propellers" are features in Saturn's A ring associat ed with moonlets tha t 



open partial gaps. They exhibit non-Keplerian motion ( ITiscareno et al.ll2010l ); 



the longitude residuals of the best-observ ed propeller, "Bleriot ," appear con- 



sistent with a sinusoid of period ~4 years. iPan fc Chiang! (120101 ) proposed that 
propeller moonlets librate in "frog resonances" with co-orbiting ring material. By 
analogy with the restricted three-body problem, they treated the co-orbital ma- 
terial as stationary in the rotating frame and neglected non-co-orbital material. 
Here we use simple numerical experiments to extend the frog model, including 
feedback due to the gap's motion, and drag associated with the Lindblad disk 
torques that cause Type I migration. Because the moonlet creates the gap, we 
expect the gap centroid to track the moonlet, but only after a time delay tdeiay, 
the time for a ring particle to travel from conjunction with the moonlet to the 
end of the gap. We find that frog librations can persist only if t delay exceeds the 
frog libration period Put,, and if damping from Lindblad torques balances driving 
from co-orbital torques. If tdeiay "C Pub, then the libration amplitude damps to 
zero. In the case of Bleriot, the frog resonance model can reproduce the observed 
libration period Put, — 4 yr. However, our simple feedback prescription suggests 
that Bleriot's tdeiay ~ O.OlPub, which is inconsistent with the observed libration 
amplitude of 260 km. We urge more accurate treatments of feedback to test the 
assumptions of our toy models. 



1. Introduction 

"Propellers" observed by the Cassini spacecraft in Saturn's A ring appear as S-like 



features superimposed on azimuthally long and radially narrow gaps ( ITiscareno et al.l 12006 
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Sremcevic et al.l l2007t iTiscareno et al.l 120081 . I201CH ) . Each propeller is believed to trace a 
moonlet several hundred meters in size which gravitationally repels nearby ring particles, cre- 
ating an underdense gap in the moonlet's immediate vicinity, as well as overdensities at radii 
just outside the gap. Because of Keplerian shear, these density perturbations propagate to- 
ward greater longitudes inside the moonlet's orbit and sm aller longitudes outside, producing 
two l obes separated by a few Hill spheres of the moonlet (ISeifi et al.ll2005l ; iLewis fc Stewart 
20091 ). The gap's azimuthal length is set by the moonlet mass a nd the time for ring parti 



cles to diffuse back in t o the gap via particle-particle interactio ns (ISpahn fc Sremcevic 
Sremcevic et~al1l2002[ ISeifi et al.llioost ILewis fc Stewart! I2OO9I ) . 
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Propellers in the outer A ring exhibit non-Keplerian motion (ITiscareno et al.ll2010l ). In 
particular, the longitude residuals of the propeller "Bleriot," observed 89 times over 4.2 years, 
show variations consistent with a sine wave of half-amplitude 260 km and period 3.68 years. 
The longitude residuals imply semimajor axis variations of order 100 m. These data constrain 
the underlying mechanism. The variations' smooth sinusoidal character suggests that this 
mechanism is not stochastic on timescales shorter than a few years. Secular interactions 
with the Saturnian moons, rings, or equatorial bulge have timescales typically much longer 
than a few years; in any case, secular interactions cannot induce semimajor axis changes. 
Finally, no other Satu rnian moon appears to occupy a mean motion resonance with Bleriot 
ITiscareno et al.lboioh . 



Pan & Chiang (2010, hereafter PC) proposed that Bleriot's non-Keplerian motion is 
caused by gravitational interactions between Bleriot's moonlet and co-orbital ring material 
outside the moonlet's gap. The interaction is long-range: the co-orbital mass is located 
thousands of Hill sphere radii away from the moonlet. The propeller moonlet and the much 
more massive co-orbital material participate in a 1:1 resonance reminiscent of tadpole and 
horseshoe orbits in the conventional restricted three-body problem. The moonlet performs 
what PC called "frog" librations within the gap. 

To re-cap the main ideas behind PC's analysis we give the following order-of-magnitude 
description of frog librations. Because our goal here is to describe the physical processes 
clearly, we neglect factors of order unity; we refer readers to § 2 of PC for an exact treatment. 
We use units where Saturn has mass Mg aturn = 1, the co-orbital mass's distance from Saturn 
is r coor b = 1, and Newton's constant G is 1. We work in the frame co- rotating with a test 
particle on a circular orbit of radius 1 and (as a consequence of our unit system) angular 
velocity 1. The moonlet of mass /i m00 n opens a gap of angular size 20 <C n in the co-orbital 
material (Figured]). The co-orbital material outside the gap has IXlcLSS /-/-ring* 

Since the co- 
orbital material closest to the ends of the gap interacts most strongly with the moonlet, we 
model the co-orbital material as two point masses, each of mass /i en d = Airin g 0/2, fixed at 
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longitudes ±0. We label the moonlet's polar coordinates (r — 1 + A, 8) where A 1 and 
\9\ < The moonlet librates about the equilibrium pointo between the co-orbital masses, 
as shown in the right-hand side of Figure [TJ 

The moonlet's motion is governed by Keplerian shear and gravitational interactions 
with the co-orbital masses. We derive scaling relations for the frog libration's period and 
radial/azimuthal aspect ratio as follows. The azimuthal length # max of one libration is of 
order the moonlet's azimuthal Keplerian drift in one libration period Put,. If the libration 
radial width is A max , then the moonlet's drift speed is of the same order, and 

#max ~ -Plib " A 

max • V J 

Over the libration cycle, gravitational interactions with the co-orbital masses change the 
moonlet's semimajor axis from 1 + A max to 1 — A max . The corresponding fractional change 
in the moonlet's specific angular momentum is ~A max . The angular momentum changes 
because the (azimuthal) gravitational forces ~±/z en d/0 2 exerted by the co-orbital masses do 
not cancel exactly; the residual is of order 



Together, Eqs. [1] and [2] imply 



» /^end $max D / s 
A max ~ —7, — ■ "lib • 



r> $ ^ Amax ^end (o\ 

These scalings match those of the exact relations in Eqs. 10 and 11 of PC, which we repeat 
here for reference: 

Note that the coorbital mass [i used by PC equals 2/i en d as we define it here. 

In a more detailed calculation with the co-orbital mass spread uniformly over the entire 
co-orbital region outside the gap, PC confirmed the scalings of Eq. [3] and found that for 
parameters characteristic of Bleriot's environment, the frog librations should indeed have a 
period of ~4 years. 

While the agreement between the period of the frog reson ance as estimate d by P C, 



and the period of Bleriot's longitude variations as observed by iTiscareno et al.l (120101 ) is 



1 This equilibrium point is a local maximum of the gravitational potential. The moonlet can librate about 
this potential maximum because the Coriolis force stabilizes its motion. 
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Fig. 1. — Schematic of frog librations. The left side shows the geometric configuration of the 
moonlet and its co-orbital material. We use units where Saturn's mass is unity and where 
the orbit radius of the co-orbital ring material (dashed circle arc) is unity. The moonlet 
of mass /imoon located at polar coordinates (r = 1 + A, 0) moves in its gap of angular size 
20 in the co-orbital material. The portions of the co-orbital mass /x r i ng that interact most 
strongly with the moonlet are those within azimuth ~0 of the ends of the moonlet's gap. 
For simplicity, we model the entire co-orbital mass as two identical point masses of mass 
/Ucnd = Urmg^/Z located at the gap ends. The dashed oval curve on the right represents a 
frog libration trajectory. 
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promising, PC's treatment is incomplete. In particular, PC treated the co-orbital masses 
as stationary (in the co- rotating frame). This is a severe simplification. Because the moon- 
let creates and maintains its gap, the co-orbital material must, to some degree, follow the 
moonlet's position. If it follows too closely, the frog libration amplitude (8) may be unob- 
servably small. On the other hand, if the gap ends are far enough away from the moonlet, 
the mass there will respond sluggishly to the moonlet's motion because of the finite time 
needed for particles to drift from conjunction with the moonlet to the gap ends. In this case, 
large-amplitude frog librations should be permitted. 

A self-consistent understanding of the moonlet's motion must account for how the non- 
Keplerian motions of the moonlet feed back into the non-Keplerian motions of the co-orbital 
masses. Here we investigate the behavior of frog orbits when the co-orbital masses move 
in response to the moonlet, and when Lindblad torques due to close encounters between 
the moonlet and ring particles are present. Our focus here is on conceptual clarity; from a 
technical standpoint, our models are, by design, crude. Many of our calculations are accurate 
to order-of-magnitude at best. In our equations we use "~" to denote a relation in which 
order-unity factors are ignored; "~" to denote a relation in which order-unity factors are 
retained but some quantities remain approximate or not precisely defined; and "=" to denote 
an exact relation. In Section 12.11 we present a simple toy model for feedback and examine 
its implications for frog librations. In Section 12.21 we add Lindblad torques. In Section [3] we 
summarize and comment on the results of our experiments. 



2. Modeling the gap motion and Lindblad torques 

Figure [2] gives an overview of the formation and maintenance of the moonlet gap. The 
moonlet sits in a sea of much smaller ring particles with mass surface density E. In close 
encounters with ring particles, the moonlet gravitationally repels the particles, clearing a 

1/3 

region whose initial radial size is of order the moonlet's Hill sphere /imoon- As repelled 
particles drift away in longitude from the moonlet because of Keplerian shear, they collide 
with other particles; over a time ideiay, they diffuse back into the gap, refilling it. Over 
this diffusion time, the perturbed ring particles drift an angular distance from their close 
encounter with the moonlet to the end of the gap. That is, 

_ 2 

^delay — 7T (5) 
"J *^gap 

1 /3 

where x gap < //moon is the gap's typical downstream radial half-width (radial distance between 
the moonlet's orbit and one edge of the gap). In this definition x gap is the gap width along 
most of the gap's (azimuthal) length, far downstream of the moonlet's location; its value 
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depends on the way ring particles diffuse back into the gap and, therefore, on the local ring 

1 /3 

viscosity. The downstream gap width x gap may differ from the gap width x ~ //moon at the 
moonlet's longitude. 

We can think of ideiay as the delay between any movement of the moonlet and the 
response of the gap ends to that movement. If the moonlet were to suddenly shift its position 
in the rotating frame, the gap ends would re-position themselves so as to be centered on the 
moonlet's new position — but only after a time tdeiay This simple picture suggests that we 
modify the model of PC so that at any given instant, the co-orbital masses are located ±0 
in longitude away from the moonlet's location a time ideiay ago. 

2.1. Experiments including gap motion only 

To explore the effects of this delayed gap motion on the moonlet, we construct a numer- 
ical toy model. As in PC, we treat the co-orbital material as two point masses each of mass 
/i en d located at the ends of the moonlet gap, and we integrate the trajectory of the moonlet 
(test particle) as it interacts with the point masses and the central unit mass. To model the 
feedback, we fix the point masses' angular separation to be 20 and their radial coordinate 
to be unity, and force them to move so that the longitude of their midpoint at time t equals 
the moonlet's longitude at time t — tdeiay To begin the calculation, we integrate from t = 
to t = tdeiay with stationary point masses, and then use this history when we turn on the 
feedback at time t = ideiay We treat /i en d, 0, and ideiay as free parameters. 

The top panel of Figure E] shows representative results from this toy model. When tdeiay 
is short compared to the nominal frog libration time (Eq. UJ), the librations damp quickly. 
When tdeiay is comparable to p^ or longer, the libration amplitude grows and the moonlet 
eventually exits the resonance. The transition between amplitude growth and amplitude 
decay occurs when P lih is a few times t delay ; in the representative example shown in Figure El 
this transition occurs when tdeiay ~ Pib/3. 

The shape of the potential the moonlet experiences also changes because of the co-orbital 
masses' motion. The change in the potential alters Put,: in Figure [3j P u b — 60 rather than 
Pub — 81.6 as predicted by Eq. HI In our numerical trials using parameter values spanning 
the ranges 10~ 7 < // cnd < 10~ 4 , 10~ 3 < < 0.5, and 10 < tdeiay < 7000, we found that 
the ratio between tdeiay and P^ is the principal deciding factor for the moonlet's qualitative 
behavior. 

We interpret the growing librations when tdeiay ^ Pib as resonant forcing. In this case, 
the gap position lags the moonlet position by a phase of order unity or, equivalently, a 
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>- -<^-Plib ^delay 
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direction of orbital motion X ^ ^ _ qO 




Fig. 2. — Schematic showing the trajectory of a ring particle (shaded circle with arrow) 
participating in gap formation. A moonlet of mass /i m00 n is embedded in a particle disk of 
surface density S. For a particle with initial radius inside the moonlet 's orbit (toward the top 
of the figure), a gravitational interaction with the moonlet decreases the particle's semimajor 
axis. Diffusion via collisions with nearby ring particles moves the particle back to its original 
semimajor axis over a time tdeiay During this time the particle drifts an angle relative to 
the moonlet because of Keplerian shear. The dashed oval represents frog librations of the 
moonlet with period P\^. 



-8 - 



longitude of order the libration amplitude. Because the co-orbital masses are slaved to the 
moonlet's actual past motion, they automatically drive the moonlet at resonance^] 

In the limit tdciay <^ P\ib, however, the moonlet's displacement from the gap center is 
always small compared to the libration amplitude that a moonlet with the same location and 
velocity would have if the co-orbital masses were stationary. The moonlet's speed relative 
to the gap is also much less than the speed it would have if the co-orbital masses were 
stationary. Continuously shrinking the moonlet's displacement from the equilibrium point 
between the co-orbital masses, and its speed relative to the equilibrium point, damps the 
libration amplitude until the moonlet lands at the fixed point of the potential. 

We find that the transition between the regimes of growing and decaying libration am- 
plitude occurs over a narrow range in tdeiay For example, the values for tdeiay corresponding 
to the two different outcomes in the top panel of Figure [3] differ by less than 10% of Pm,. 

We performed some supplementary experiments. First we forced the co-orbital masses 
to move sinusoidally in longitude with prescribed amplitude and frequency. Setting the 
co-orbital masses' oscillation period equal to P^ as given by Eq. @] produces beats in the 
moonlet's motion: the motion of the co-orbital masses changes the shape of the potential 
the moonlet sees, shifting the libration period enough for the driving to be just off resonance 
(Figure H]). If we decrease the co-orbital masses' oscillation amplitude while keeping their 
oscillation period fixed, the potential experienced by the moonlet changes less. Then the 
driving is closer to resonance and the beat frequency decreases. Compare these results with 
those of our original experiments, where the co-orbital masses were slaved to the moonlet's 
actual past motion and the driving was therefore exactly on resonance. 

In experiments aimed at softening the transition between the regimes of growing and 
decaying amplitude, we performed simulations where we artificially weakened the feedback 
strength. To do this we forced the co-orbital masses to follow the moonlet's past motion, as 
before, but scaled the amplitude of the co-orbital masses' motions relative to the moonlet's 
by a constant factor. We observed no qualitative change in the moonlet's behavior. The 
value of tdeiay at the transition changed by a factor of order unity. 

To summarize our findings so far: the libration amplitude damps to zero if ideiay is too 
short and grows without bound if ideiay is comparable to or longer than p^, with an abrupt 
transition between the two regimes. However, longitude residuals like those observed might 



2 Since the co-orbital masses follow the moonlet's past rather than its current motion, the driving frequency 
may, in principle, differ slightly from resonance, particularly for idciay 3> -Piib- However, the driving frequency 
always moves toward resonance with time. 
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be produced in the latter regime if some other mechanism damps the motion enough to limit 
the libration amplitude's growth. We explore one possible damping mechanism in §2.21 



The requirement that tdeiay ^ -Pnb/3 that we find in our simulations (the factor of 3 
is empirical) has a simple consequence for the gap's radial width: using Eq. @] and /i en( j — 



X • 2x gap ■ 0, we find that 



7T 3 / 2 7T 



-gap 

which together with Eq. [5] implies that 



or 



< 4 m 



x g&p <- 24 Sr coorb 

'"coorb -^Saturn 



'"coorb 



^gap ~ ^ 4Q g / cm 2 J y iMx 1Q 10 cm ^ 

where for t he last line we have re s tored the units and used numbers inspired by Bleriot and 



the A ring (ITiscareno et al.ll2010l ; IColwell et all 120091 ). 



At face value, Eq. [S] suggests that for Bleriot's moonlet to avoid libration damping, its 
gap should be radially narrow, of order several meters across, over a significant portion of its 
azimuthal extent. The required narrowness seems implausible, as it is compara ble in size to 



the m eter-sized particles dominating the mass in the A ring (see Table 15.1 of ICuzzi et al. 



20091 ). Thus, Eq. [8] argues against frog librations as being the correct explanation for the 
observed non-Keplerian motions of propellers. We stress, however, that our equations are 
only as accurate as the toy model for tdeiay from which they derive. A more careful analysis 
of the gap's reactions to the moonlet (e.g., using numerical simulations of the entire gap) 
would provide an important test of our estimates here. 



2.2. Experiments including gap motion and Lindblad torques 

The results of §2.11 indicate that stable frog librations can persist only if tdeiay ^ -Piib 
and if some other mechanism damps the librations enough to stop their runaway growth. 
One damping mechanism is disk torques due to Lindblad resonances — the same torques that 
open the gap and are responsible for Type I migration. Because Lindblad torques repel the 
moonlet away from the gap edges, forcing the moonlet's semimajor axis back toward the gap 
center, we expect such torques to decrease the frog libration amplitude. 

To provide a conceptual basis for the behavior we expect in our numerical experiments, 
we give the following order-of-magnitude estimates of the Lindblad torques and associated 
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Fig. 3. — Sample moonlet trajectories showing the effects of feedback from co-orbital mass 
motion and damping from Lindblad torques. The parameter values /i en d = 10 -4 ' 5 and <ft — 0.4 
are used throughout, so that the nominal libration period according to Eq. H]is Pub = 81.6. 
Including co-orbital mass motion changes the potential the moonlet sees; in all the sample 
trajectories shown here, the libration period in the altered potential is actually ~ 60. 
The top panel shows the effects of co-orbital mass motion alone with different values of ideiay, 
the time by which the co-orbital mass motion is delayed relative to the moonlet's motion (i.e., 
the gap re-adjustment time). Values of t delay greater than some threshold lead to increasing 
libration amplitudes (solid curve, tdeiay = 30). Values of tdeiay even slightly smaller than 
the threshold produce decreasing amplitudes (dotted curve, tdeiay = 25). The bottom panel 
includes the effects of both gap movement and Lindblad disk torques. All three trajectories 
have tdeiay = 30 and £ = 0.01, but the moonlet masses are /x m00 n = {1-48, 7.41, 37.1} x 10~ 7 , 
respectively, for the light solid, heavy solid, and dotted curves, and the gap widths are x = 6 
moonlet Hill radii . The coefficient of 6 is giv en approximately by the edges of the propeller 
gaps computed by lSpahn fc Sremcevid (120001 ) . The trajectories shown here are representative 
of our numerical experiments and are meant to illustrate our findings; the parameter values 
used do not correspond to the propeller Bleriot. 
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Fig. 4. — Sample moonlet trajectories from experiments where the co-orbital masses move 
sinusoidally with prescribed amplitude and period. For both panels, the prescribed period 
is 81.6 (equal to the period that the moonlet would have if the co-orbital masses were 
stationary). The top panel shows the moonlet's motion when the co-orbital masses move 
with amplitude 0.005 and the bottom panel shows results when that amplitude is 0.05. The 
stronger driving in the bottom panel trajectory leads to both a larger beat amplitude and 
a higher beat frequency. The latter implies that the co-orbital masses' motion changes the 
potential seen by the moonlet and the natural frequency of oscillations in that potential. 
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damping timescale, dropping fac tors of order u nity. According to the standard impulse 
approximat ion (see, for example, iDermottl ( 11984 ) for a derivation using the impulse approx- 



Goldreich &: Tremainel ( 1982 ) for a derivation using Lindblad resonances), the 



imation or 

moonlet's azimuthal acceleration due to its "one-sided" interaction with ring material on 
one side of its gap is of order 

j-i MmoonS / n \ 

0,1-sided ~ n ■ (yj 



X" 



This is consistent with the scaling A J oc /i moon /a; 5 for the one-sided impuls e of angular 
mome ntum per unit disk mass received by the moonlet (see, e.g., Eq. 32 of iCrida et al. 
( 120101 ): multiply their equation by the disk mass ~Sa; 2 that passes conjunction with the 
moonlet per unit time — this is analogous to the integral over x that they describe with 
their Eq. 36 — and divide by /z m00 n to get an acceleration rather than a force). Note that 
we use x here rather than x gap because we are concerned with ring material that just passes 
conjunction with the moonlet, i.e., ring material near the moonlet's longitude. 

The moonlet interacts simultaneously with ring material radially interior and exterior 
to it. If we assume the ring surface density is uniform, then the two one-sided torques cancel 
up to any asymmetry in the moonlet's radial distances to the gap edges. Then the moonlet's 
net acceleration is 



A 

• -T0,l-sided ~ 

X 



XA 



X 1 



(10) 



The minus sign enters because the moonlet is repelled more by whichever gap edge is closer 
to it. The damping timescale is 



''damp 



e_ 



X 



(11) 



where we have used | A | ~ \6\ due to Keplerian shear (see Eq. [T]) . 



We combine Lindblad damping with co-orbital feedback by adding the term Fg (as 
given by Eq. [TUJ with an assumed coefficient of unity) to the moonlet's azimuthal equation 
of motion in the numerical toy model of Section 12.11 We expect that damping will balance 
the (positive) feedback when tdamp is of order the libration growth timescale. In the example 
shown in the top panel of Figure [31 the timescale to amplify librations is about t a mpiif y = 
6 libration periods of the undriven system (driving alters the potential and shortens the 
libration period). To check quickly that we understand our simulations in the bottom panel 
of that figure, we make a simple order-of-magnitude estimate. We set tdamp = ^amplify = 6Piib 
and substitute Eqs. HJ [TTJ and the parameters of the example in Figure |3j /i en d = 10~ 4 5 , 
x = 6(/i m00 n/3) 1 ^ 3 or six Hill radii of the moonlet, X = 0.01, and <p — 0-4. In this example, 
damping balances growth when /i m00 n ~ (3 5 ^ 2 vr 3 /6 9 )-0 9 / 2 S 3 /i^ i / ' 2 ~ 4x 10~ 6 . Our simulations 
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show that a somewhat smaller /i m00 n ~ 7.5 x 10 7 is needed (bottom panel of Figure [3]), but 
the agreement is reasonable given the order-of-magnitude nature of our arguments. 

Since our Eqs. [TU1 1 hrough ITT1 ne gleet all order-unity coefficients — strictly speaking, they 
are scaling relations — we may apply them to real systems only with the understanding that 
numerical results cannot be more than broadly suggestive. If we use x ~ 6 moonlet Hill radii 
as before, then Eqs. HI HU and the equality in Eq. [7J (x gap /S ~ 24/n 2 ) give 

^damp q g / faoon \ ( 4> \ (12) 



P, ih ' V 3 x 10- 15 J V 0.004 



where the numerical values are inspired by Bleriot ( ITiscareno et al.ll2010l ). This numerical 
result suggests that positive co-orbital feedback and damping Lindblad torques might plau- 
sibly balance for propellers in the outer A ring. Unfortunately, the result relies on Eq. [Jj 
which as we have seen previously demands an unrealistically narrow gap. 



3. Summary and Discussion 

As a step toward a self-consistent model for propellers' non-Keplerian motions, we added 
to the frog libration model of PC new terms representing co-orbital mass motion and Lind- 
blad torques. In numerical experiments with this 'extended' frog model, we found that 
allowing the co-orbital mass at the ends of the gap to move in response to the moonlet 's 
non-Keplerian motions can either damp the frog librations completely or drive the librations 
resonantly. In other words, in our simple model for co-orbital feedback, the feedback is either 
strongly negative or strongly positive. We found further that the strong positive feedback 
could be limited, and the motion stabilized, by Lindblad torques. 

Our numerical experiments emphasize simplicity over realism. Although they clarify 
the conditions necessary for stable frog librations, they offer no explanation for why these 
conditions should be met. That is, we leave unanswered the question of why tdeiay, -Pub, 
/-tmoon, and /i r i n g of real propeller systems should occur in the right proportions for stable 
frog librations to exist. 

Still, our extended frog model makes observationally testable predictions: if propeller 
moonlets are performing frog librations, we expect that 1) the properties of Bleriot 's gap 
(its azimuthal and radial dimensions, and the co-orbital surface density outside the gap) are 
such that its frog libration period Pub — 4 yr; 2) Bleriot 's and other propellers' longitude 
residuals will continue to vary sinusoidally in time; 3) at any given time, propeller positions 
will be typically offset from the centers of their gaps by azimuthal distances of order the 
observed rms longitude residuals; and 4) the properties of propeller gaps are such that gap 
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drift (equivalently, gap closing) timescales tdeiay (Eq. |5]) are comparable to or longer than 
frog libration periods P\^. Of these four predictions, the first three are retained from PC's 
original frog model, while the fourth is new from the extended frog model. 

Unfortunately, we cannot seem to satisfy both predictions 1 and 4 for Bleriot. If we take 
the azimuthal length of Bleriot's gap to be ~ 0.004 and its radial width x gap to be a few km 
(a few moonlet Hill radii), we can indeed reproduce P^ ~ 4 yr, in accord with the Cassini 
observations (PC, see the discussion following their Eq. 12). But a gap of such dimensions 
would close in a time tdeiay ~ O.OlPub (Eq. [5]) and would lead to a libration amplitude much 
smaller than the 260-km amplitude that is observed. For the gap drift time tdeiay to be longer 
than the frog libration time Put,, the radial width of the gap may need to be unrealistically 
small - - comparable to the ring particle size. This conclusion is tentative because our 
treatment of feedback is primitive and possibly oversimplified. We state this shortcoming of 
our model in the hope that more accurate studies — e.g., numerical simulations designed to 
explore long-range interactions between the moonlet and the entire gap — may either confirm 
that the frog resonance is not responsible for the observed non-Keplerian motions, or reveal 
that feedback works in a way different from what we have imagined in this paper. 

We thank an anonymous referee for a detailed reading of this work. 
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